Effects of Precipitation on Propagation at 
0.63, 3.5, and 10.6 Microns 

By T. S. CHU and D. C. HOGG 

The attenuation and scattering of laser beams by rain, fog, and snow 
have been calculated and measured at 0.63, 3.5, and 10.6fi. Attenuation 
of the infrared wavelengths by light fog is up to one order of magnitude 
less than at 0.63n. But for dense fog, calculation shows that the attenuation 
at 10. 6 n can exceed 40 dB per km. It is found that attenuation by rain 
can be calculated to good accuracy from average path rain rates provided 
that forward scattering is taken into account; this scattering reduces the 
attenuation. 

Measurements of propagation through precipitation over a 2.6 km path 
are discussed in detail and are found to be consistent with predictions. 
The wavelength dependence of attenuation is found to vary from one fog 
to another because of different drop-size distributions. Attenuation of 0.63\i 
in rain showers is about 20 per cent less than at 8.5y. and is, of course, 
much less than the attenuation in fog. But even for extremely heavy rain 
showers, the 0.63n attenuation never exceeded 20 dB per km, which is 
less than the attenuation of millimeter waves under such conditions. Both 
the attenuation and foreward scattering properties of snow appear to be 
between those of fog and rain. 

I. INTRODUCTION 

It is well known that optical waves suffer much attenuation in 
propagation through precipitation in the atmosphere. Penetration of 
light through a dense fog is much more difficult than through a heavy 
shower. One would also expect similar phenomena for near infrared 
waves. However, quantitative data are necessary in order to evaluate 
the potential use of lasers for short range communication through 
the atmosphere. Attenuation by precipitation is the dominant factor 
in determining the feasibility of open-air laser communication sys- 
tems. The interesting fluctuations of laser beams caused by turbulence 
in the clear atmosphere are not discussed here. 
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The gas lasers which provide significant amounts of power within 
the infrared windows of the clear atmosphere are 3.5/x and 10.6/x; 
3.39/x is not suitable because it is attenuated by the small amount of 
methane in the air. 1 The wavelengths of visible gas lasers with suitable 
output power are 0.49^ and 0.63/1. Since the attenuation of optical 
waves in a fog is largely the result of scattering by numerous small 
drops, it is not likely that transmission of 0.49/x and 0.63^ through a 
fog would differ appreciably. The literature reveals that the attenua- 
tion of blue light (~0.49/i) through a light fog is somewhat greater 
than that of red light (~0.63/x), while their attenuations in a dense 
fog differ little from one another. 2 In this study, 0.63/a is the visible 
wavelength used. 

Some computed results of attenuation by precipitation at the wave- 
lengths of 0.63/a, 3.5/*, and 10.6/i are given first. The effect of forward 
scattering is examined using a Gaussian model for the scattering pat- 
tern of a rain drop and a single scattering model is used to describe 
the broadening of a laser beam by precipitation. Then the results of 
measurements of attenuation by precipitation over a 2.6 km path at 
these three wavelengths are described. Some of the measured data 
have been discussed briefly in previous notes. 3, 4 

The measured attenuation by rain is compared with that calculated 
from the rainfall measured on the path. Optical attenuation by fog 
is compared with total liquid water content derived from path-loss 
measurements taken on a radar operating at a wavelength of 3 mm. 
Similar quantitative correlations have not been made for attenuation 
by snow because of the lack of an independent physical measure- 
ment for this form of precipitation. 

II. ATTENUATION CALCULATIONS 

2.1 Method oj Computation 

The total power removed from a laser beam by a precipitation par- 
ticle, say a spherical drop of radius a, is given by 

E = SC(a, X) (1) 

where S is the magnitude of the incident Poynting vector, and C(a, A) 
is the total or extinction cross section of the drop. The change in S 
upon passing through a layer of precipitation Az is 

AS = -S Az [ n(a)C(a, X) da (2) 

Jo 
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where n(a)da is the number of drops per unit volume within the range 
(a, a+da). Integrating equation (2) yields 

8 = S n e- J ' ad! (3) 

where 

a = f n(a)C(a, X) da. (4) 

Jo 

The total cross section is often normalized with respect to the geo- 
metrical cross section 

Q(a, X) = S2A (5) 

7ra 

Since the energy removed is partly scattered and partly absorbed by 
the particle, the total cross section is divided into scattering and 
absorption cross sections: 

Q.xt = Qscat + Qub, • (6) 

The Q's are sometimes called efficiency factors. Notice that part of 
the energy represented by the scattering cross section is scattered 
forward ; this leads to a decrease in attenuation as discussed in Section 
III. 

The attention coefficient in equation (4) may be calculated if the 
particle-size distribution and the cross sections are given. Equation 
(4) gives the extinction, scattering, or absorption coefficients when the 
appropriate cross section is substituted.* We take the shape of the 
water drops in fog and rain to be spherical. With this assumption, it 
is not a prohibitive task to use Mie's rigorous expressions to calculate 
the extinction and absorption cross sections of the drops. However, 
very simple approximate formulas are quite applicable here in view 
of the uncertain drop-size distributions and the doubtful sphericity. 
Using a combination of geometrical optics and Huygen's principle, 
Van de Hulst 5 obtained the following formulas for the extinction and 
absorption cross section: 

n o i -Pt-ng/' cos/^ . . 

Q ext = 2 - 4e p l 1 sin (p - 0) 

, -u tanfl/cOS /3\" , __. , ./COS/SY" „„ ,„ „ 

_ A( , P«u^ J cog(p _ 2j8 j + 4 ^ Ml cog2 ^ ^ 



* Extinction is often used in optical terminology for attenuation; however, it 
is misleading. It is more exact to designate the extinction coefficient as the at- 
tenuation coefficient without forward scattering corrections. 
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• Qab, = 1 + 

where 



2xri 



+ 



- 1 



8zV 



(7b) 



Aura , , 

P = — (n - 1) , 



tan /3 = 



n - 1 ' 



JE = 



2wa 



and n and ra' are the real and imaginary parts of the refractive index, 
respectively. The above two expressions are derived under the assump- 
tions (n — 1) « 1, n' « 1, 2ira/\ » 1. Numerical comparison indicates 5 
that the essential features of the cross section curves using equation (7) 
remain correct for | n - 1 - in! \ < 1 and 2ira/\ > 1. When p « 1, 
equation (7) reduces to 



Q„, = ip tan + \p\l ~ tan 2 j8) 
Q„ ba = ^P tan p - p 2 tan 2 j8. 



(8a) 
(8b) 



2.2 Refractive Indices 

Figure 1 shows the complex refractive index (n — in') of liquid 
water in the visible and near infrared region. For the wavelengths 
below 18/a the data have been taken from the tabulation by Centeno 
except for those of n' in the 2.8/* band where Centeno's values are too 
low and have been replaced by the Bayly results. 7 No serious changes 




0.1 0.2 0.4 0.6 1.0 2 4 6 10 20 40 60 100 200 

WAVELENGTH IN MICRONS 



Fig. 1 — Complex refractive index (n — in') of liquid water. 
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at other wavelengths are reported in more recent literature. In par- 
ticular, for the three wavelengths we investigated, the complex re- 
fractive indices of liquid water are well agreed upon : 



X = 0.63m, n = 1.33, 



X = 3.5/*, 



n = 1.42, 



ri ~ 10"" 
n' = 0.013 



X = 10.6m, n = 1.18, n' = 0.08. 

The curves beyond 18/x in Figure 1 have been obtained from Kislov- 
skii. 8 Most measurements of the refractive index of water have been 
conducted at 20°C, but it is known that as the temperature of the 
water is increased, the position of maximum absorption is shifted to 
slightly shorter wavelengths. 0, 10 However, detailed information on 
temperature dependence seems to be lacking. Much less data have 
been published for the complex refractive index of ice than for water. 
Figure 2, ice at — 10°C, is reproduced from Kislovskii for the sake of 
completeness. 

2.3 Cross Sections and Attenuation Coefficients for Mo?iodisperse Drop- 
Size Distribution 

Using equations (7) and the above values of the complex refractive 
indices, extinction and absorption cross sections have been computed 
as a function of radius of water spheres for the three wavelengths 
we mentioned. The data are shown in Figure 3. When the radius of 
the water sphere becomes large, say 1 mm (that is, about the size of 
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Fig. 2 — Complex rofrartive index (n — in') of ire. 
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6 8 10 12 14 16 

RADIUS OF WATER SPHERES IN MICRONS 



Fig. 3 — Approximate extinction and absorption cross sections. 

a rain drop), the extinction cross sections approach two, while the 
absorptive cross sections at 3.5/x and 10.6/x approach unity. The ab- 
sorption cross section of 0.63/* does not appear in Figure 3 because 
pure water is practically lossless at visible wavelengths. The extinc- 
tion cross section for 3.5/x remains substantially less than that of 0.63/x 
only when the radius of water sphere is smaller than lju.; also, in that 
range, the extinction cross section of 10.6/* almost entirely results from 
absorption. The longer wavelength and the lower refractive index push 
the first maximum of the 10.6/t extinction cross section to an "a" of 
about 19ti. 

Although fogs and rains never have drops of a single radius, it is 
instructive to calculate the attenuation coefficients for this ideal case. 
If we substitute 



w 



n(a) = t — s <5(er) 
s-n-a 



00 



into equation (4), where w is the liquid water density and 8(a) is the 
Dirac delta function, then the extinction coefficient of liquid water 
becomes 



3 Qext nepers /mg = g ^ Q^ dB / mg 

4 a km / m 8 a km/ m 3 



(10) 
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where the radius of water sphere "a" is in microns. The extinction 
and absorption coefficients obtained using equation (10) are plotted 
in Figure 4. 

The vast difference in attenuation by rain and fog is very clear 
from equation (10). Although the liquid water content of a typical 
heavy shower, say 1000 mg/m 3 , is ten times that of a typical dense 
fog, say 100 mg/m 3 , the fog drop radius is about one thousandth that 
of a rain drop. This is the primary reason why attenuation by rain 
is two orders of magnitude less than that of fog at optical wave- 
lengths. For rains with drop radii about 1mm, the extinction coef- 
ficient is only 0.0065 (dB/km)/(mg/m 3 ) (see Figure 4), thus a 
fairly heavy shower with a liquid water content of 1 g/m 3 (rain rate 
about 25 mm per hour) has an extinction of about 6.5 dB/km for 
both optical and near-infrared wavelengths. 

At least half of the attenuation resulting from rain at 3.5ju and 
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Fig. 4 — Extinction and absorption coefficients of precipitation with mono- 
disperse drop-size distribution. 
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10.6/*. is absorptive loss while at 0.63/a it is caused almost entirely by 
scattering, as seen from Figure 4. 

Attenuation by fog is contributed largely by water drops less than 
a few microns in radius, therefore much of the 10.6/x attenuation is 
absorptive loss, as implied by equation (8) and Figure 4. The extinc- 
tion coefficient of fog is expected to be substantially less (up to one 
order of magnitude) at 10.6/* than at 0.63/x where scattering predomi- 
nates. However, heavy fogs with liquid water contents exceeding 
O.lgm/m 3 will attenuate 10.6/* more than 40 dB/km. 

Attenuation by fog at 3.5/a, largely caused by scattering, remains 
significantly less than that of 0.63/t only if the fog drop radii are 
smaller than 1/x, as is also true for haze. Fogs which consist of drops 
with an average radius of 5/i are expected to attenuate about the 
same amount at optical and shorter near-infrared wavelengths. 

2.4 Calculations for Realistic Drop Size Distributions 

Experimental difficulty and natural variations have caused gross 
inconsistency among the published data on measured fog drop-size 
distributions. Most of the measuring instruments fail to respond to 
droplets less than one micron in diameter, yet there are indications 
that the number density continues to increase as the diameter of 
fog drops decreases toward one micron. The relatively consistent 
information seems simply to be that the number density decreases 
for the larger drop radii. 

For example, Eldridge's measurement shows that the number den- 
sity decreases approximately with the inverse square of the drop 
diameter. 11 In searching for a realistic drop-size distribution for fog, 
the following model 12 for the number density, appeared to be plausible: 

Here 

A = (*)" -r^—Z and B-S. 

N is the number of drops per unit volume (say cm 3 ), Om the radius 
of the drops with the maximum number density, and a and y, param- 
eters which can be adjusted to fit an observed fog particle size dis- 
tribution. The liquid water density appropriate to the above model 
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is easily obtained by integrating equation (11). 



IT = / n^ira 
Jn 



da = 



■W 



a + 4 



+ 1 



al,N. 



(12) 



In Figure 5 four plots of equation (11) are given for several indi- 
cated combinations of parameters. Three of these exhibit a slope on 
the larger drop side of the maximum (a > a m ) similar to that of 
Eldridge's measured distribution. The other combination of param- 
eters (a = 2, y = 3) shown in Figure 5 corresponds to the asymptotic 
distribution of coagulating particles 13 and has a sharp decreasing rate 
on both sides of the maximum number density. 

Using these model distributions, extinction coefficients have been 
plotted versus the radius at maximum number density, a m , in Figure 
6. These curves merge at somewhat smaller values of a m than do the 
corresponding curves for constant drop radius in Figure 4 because 
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Fig. 5 — Drop size distribution of fog models. n(r) = Afexp[(— a/y)r y ] A = 
(N/a m ) (a/7) (0+1)/7 y/Tla+l)/y1, and a,„ = radius of drops with maximum 
number density. 
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the drops larger than a m contribute more to the extinction than do 
smaller ones.* The similarity between Figures 6(a) and 6(c) is 
attributed to the similar average slopes on the larger drop side of 
the corresponding distributions in Figure 5. 

It is interesting to compare the above model with the Junge dis- 
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Fig. 6 — Extinction of coefficients for various fog models. 

tribution 14 which states that n{a) = ca^ where 2 < p < 4. 15 - 1G This 
simple relation points out the importance of decreasing number den- 
sity for the larger drops; however, more parameters, such as those 
offered in equation (11), are needed to describe both ends of the 

* Although more extinction is obtained jor the same amount of liquid water 
in the form of smaller drops than in the form of larger drops, the extinction cross 
section of a larger drop is, of course, greater than that of a smaller drop. 
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distribution. Since the extinction cross section is of the order of the 
geometrical cross section Tra 2 , the decreasing rate of larger drops can 
be approximated by the inverse square power (p <•* 2) only in a 
limited range beyond which the number density must decrease much 
more rapidly to insure convergence of the integral in equation (4). 
Raindrop size distributions measured at ground level seem to be 
fairly consistent among the published data. The Laws and Parsons 17 
distributions which have been virtually confirmed by later workers 
are illustrated in Figure 7 and are readily available in tabulations. 18 
Substituting these well-known distributions into equation (4) yields 
the extinction coefficient versus rain rate shown as the upper solid 
curve in Figure 8. Notice that the total extinction is the same for 
the visible and near-infrared wavelengths. Half of the extinction 
coefficient at 3.5/x. or 10.6/* is contributed by the absorptive loss as 
illustrated in Figure 8. 

III. FORWARD SCATTERING 

The energy intercepted by precipitation particles is partially scat- 
tered forward, toward the receiver, and this reduces the attenuation. 
Here is an approximate calculation for predicting the reduction of 
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Fig. 7 — Raindrop size distributions by Laws and Parsons. 
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attenuation at 0.63/i and 3.5/* by forward scattering during rain, 
and a single-scattering model for beam broadening. 

3.1 Approximate Correction Factor jor Attenuation by Rain 

The forward scattering pattern of a rain drop is assumed to be a 
Gaussian beam with an effective aperture area equal to the geo- 
metrical cross section of the water sphere. 19 Using the geometrical 
parameters illustrated in Figure 9, the scattering pattern of a rain- 
drop of radius a situated at (p, <p, Z) is taken to be 

_2_ -2RV».' 

2 e 

irw. 



where 



«!si0l[i + x 



( HL - Z) 
\ twI. 



and the effective aperture area is 2iriv 2 ol = wa 2 . The beam from the 
transmitter, T, is taken to be 



■KW 




E 20 

2 ,6 



Q 8 



CALCULATED ATTENUATION 

WITH 3.5/1 FORWARD 
SCATTERING CORRECTION 



CALCULATED 0.63. 3.5 &, 10.6/1 
ATTENUATION OF THE DIRECT -j 
BEAM— TOTAL EXTINCTION 
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Fig. 8 — Rain attenuation coefficient for a 2.6 km path. 
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PRECIPITATION (p,0,Z) 

PARTICLE --^ l p A 
"P 



Fig. 9 — Geometry for forward scattering. 



where 



w = IV 



* + <£f] 



and w oT is the minimum beam radius, located at the Uniphase trans- 
mitting aperture. Then the attenuation in a path of length L is: 



aL{\ - 0) = f a 

Jo 



dZ 



0.,t 2 






TTWr 






Jo tW 


' 2 -2RVw 

zC 

TV). 


''2irp dp 



(13) 



a is the extinction coefficient (the attenuation coefficient without 
forward scattering correction), |9 the reduction factor resulting from 
forward scattering, Q BCR /Q ext the fraction of the intercepted energy 
available for forward scattering, 2/ww% the maximum flux of the 
unattenuated transmitted beam at the receiver 



k = wiA 



i + 



XL 



■kWotI _ 



and p A /(L — Z) = 9 A the acceptance angle of the receiver. It is easily 
shown from the geometry of Figure 9 that R/p = L cos 6/Z; thus 
when 6 is small, as it is in practice, R = pL/Z. If the acceptance angle 
of the receiver is large compared with the scattering beamwidth of 
the rain drop, p A can be replaced by <» . 

The above approximations simplify evaluation of the inner integral 
in (13) and lead to the following expression for the forward scattering 



736 THE BELL SYSTEM TECHNICAL JOURNAL, MAY-JUNE 1968 

correction factor: 



12- [' 

" L Q FXl J n 



w r 



w\ + </(§ 



dZ. (14) 



Using equation (14), /3 versus raindrop radius has been computed for 
L = 2.6 km, X = 3.5m, ™<>t = 0.55 cm, and X = 0.63/t, w oT = 0.25 cm. 
Table I gives the data. The two values of w oT are the clear weather 
beamwidths in our experiments. Upon numerical integration over the 
raindrop size distribution one obtains the dashed curves in Figure 8. 
The above approximate calculation is essentially a crude multiple 
scattering model for very narrow beams because equation (13) implies 
that the forward scattering energy effectively rejoins the original 
transmitting beam. 

3.2 A Single Scattering Model for Beam Broadening 

Measurements discussed later show that various degrees of beam 
broadening occur when a 0.63/x laser beam propagates through precip- 
itation. Beam broadening resulting from rain is more pronounced 
than that of snow and fog. Let us consider a single scattering model 
describing the effect of forward scattering on the attenuation, and 
beam broadening of a narrow laser beam. We are interested in the 
power distribution at the receiving plane. 

Using cylindrical coordinates with origin at T as illustrated in 

Table I — Forward Scattering Correction Factor /3 for a 
2.6 km Propagation Path 



Raindrop 


X = 0.63>i 


X = 3.5/1 


Radius (cm) 


WoT = 0.25cm 


v>„t = . 55cm 


0.025 


0.106 


0.025 


0.050 


0.202 


0.049 


0.075 


0.289 


0.071 


0.100 


0.366 


0.093 


0.125 


0.434 


0.113 


0.150 


0.495 


0.133 


0.175 


0.548 


0.152 


0.200 


0.595 


0.169 


0.225 


0.636 


0.186 


0.250 


0.672 


0.202 


0.275 


0.703 


0.217 


0.300 


0.730 


0.232 



w oT is the beam radius at which the field amplitude has fallen to 1/e 
of its maximum value in the collimated transmitting aperture. 
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Figure 9, the intersection of TP with the receiving plane is (pL/Z, 
0, L) whereas an arbitrary point in the receiving plane is designated 
by (r, <p', L). The Cartesian coordinates of these two points in the 
receiving plane are (pL/Z cos <f>, pL/Z sin <f>, L) and (r cos <£', r sin 
<j>', L) respectively. Since 6 is always very small for narrow laser 
beams, the scattering pattern of a raindrop will show negligible obliq- 
uity at the receiving plane and may be written as 



Q 2 / 2 \( P L 
S = — 2 exp < 2 I V 



~Z3 I I y cos# — r cos<£' 



)' 



irw. 



exp —2 



(#)' 



+ r 3 



to, 



+ 



+ \-y sin <j> — r sin <f>' 

I 



(15) 



2 ejjf r cos (0 - 0') 



10. 



Then the scattering by a uniform rain of constant drop radius a over 
a path of length L is : 



P. - Fa$*d* -±- f f -—^-^^Spdpdxp 

•> a Qext 2 ./„ Jo XW) 



JKWt 



(16) 



The definitions of the parameters in equations 15 and 16 are the 
same as in the preceding section. The integration over <p is recognized 
as a Bessel function of imaginary argument 



f 



exp 



4 ^- r cos (0 - 0') 



to! 



</<*> = w.(-<|gr)- (17) 



Using the well known identity 



j J (t (at) exp (-p't 2 )t dt = g-s exp 



q 

V 



(18) 



where a is an arbitrary complex number, one can perform the p-inte- 
gration in equation 16: 






J A -i 



. 4pLr\ 



p \w 2 + ¥w 



-)] 



p dp 



ztfj exp 

*{ «^V_\ fl /4£r\7 w a Z 2 w*. \1 

= 2 fcr+^zv cxp U tesy Uv + »w_ ■ (19) 
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Finally the Z-integration can be reduced to the following form for 
numerical computation: 






«'-, 



2 I 

w. + 



exp 



w 



2V 



w. 



+ (§ 



w 



dZ. (20) 



Both a and w g in the above equation are functions of raindrop radius. 
Upon numerical integration of P s over the raindrop size distribution 
one obtains the scattering pattern in the receiving plane. 

Using L = 2.6 km, A = 0.63/*, w oT = 0.25 cm (equivalent to our 
experimental clear weather half-power beamwidth of twenty seconds) 
and the Laws-Parsons distribution, scattering patterns for a light 
rain 2.5 mm per hour and a heavy rain 50 mm per hour have been 
computed; these are shown in Figures 10(a) and 10(b), respec- 
tively. The direct beam pattern, the total power pattern, and the 
normalized total power pattern are also shown there. It has been 
tacitly assumed in the computations that the scattering and the direct 
beam are attenuated by the same factor, jr - *, and that the total 
power pattern is the arithmetic sum of the direct beam and the scat- 
tering. For light rain, (Figure 10a), the scattering modifies only 
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Fig. 10 — Calculated scattering pattern and beam broadening, (a) Light rain, 
2.5 mm per hour, (b) Heavy rain, 50 mm per hour. 
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slightly the clear weather or "direct" pattern. However, in heavy rain, 
(Figure 10b), the scattering becomes predominant. 

One might expect more beam broadening the heavier the rain; 
however, since the average drop size increases with rain rate, the 
scattering pattern of a heavy rain is narrower than that of a light 
rain. This effect tends to offset somewhat the beam broadening by 
heavy rain. On the other hand, side lobes of rain drop scattering 
which have been neglected in the Gaussian approximation and mul- 
tiple scattering effects which are inherently absent in the single 
scattering model will tend to increase the beam broadening. Also 
notice that beam broadening will be quite sensitive to inhomogeneity 
of the rain. 

IV. EXPERIMENTAL EQUIPMENT 

Figure 11 is a block diagram of the transmitting and receiving 
equipment for 0.63/x and 3.5/*. Care was taken to stabilize the trans- 
mitters both thermally and mechanically. The equipment for 10.6/x 
has been described by Wilson and Penzias. 4 Figure 12 gives the 
profile of the 2.6 km propagation path at Holmdel, N. J. Table II 
lists the parameters of the transmitting and receiving systems. 

4.1 0.63n Equipment 

The 0.63/1 transmitter is a 50 mW multimode He-Ne laser; the 
cavity consists of two lOm-radius of curvature multilayer dielectric 
mirrors separated by 1.25m. The output is chopped mechanically at 
330 cycles per second and passed through a refracting telescope which 
has an input lens of 5 cm focal length and an output lens of 50 cm 
focal length; the transmitting beam has a half power beamwidth of 
about 20 seconds in clear weather. 

The received signal passes through attenuators, a 3A bandwidth 
filter, a double-convex lens with an iris in the focal plane, and a 
polarizer to a photomultiplier. The receiver acceptance angle, about 
10 minutes, is determined by the iris opening. Using a reference 
signal transmitted by telephone line, the synchronous detector gives 
a signal to noise ratio of 65 dB on a clear day; the limit of this 
sensitivity is set by the scattering of sunlight. 

4.2 S.on Equipment 

The 3.5/x transmitter is a 5 mW multimode He-Xe laser; the 1.1m- 
long cavity consists of a plane reflecting mirror and a calcite polarizer 
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Fig. 12 — Profile of transmission path. 

which acts as an output mirror. The output passes through the 330 
Hz chopper to a bisected nearfield Cassegrainian telescope, a com- 
bination of two bisected confocal paraboloids. The main reflector is 
one half of a 38.1 cm diameter spin-cast paraboloid of 30.48 cm focal 
length while the subreflector has a diameter of 1.27 cm and a focal 
length of 1.02 cm. 

Another bisected near-field telescope collects the signal for the 
receiver. The concentrated beam from the subreflector is amplified by 
two He-Xe tubes, each of about 18 dB gain. The good impedance 
match of the telescope plays an important role in successful applica- 
tion of these high-gain laser amplifiers. After the two amplifier tubes, 



Table II — Parameters of the Transmitting 
and Receiving Systems 



Wavelength (ft) 


0.63 


3.5 


10.6 


Laser oscillator 


He-Ne 


Ile-Xe 


C0 2 


Transmitter power (mW) 


50 


5 


500 


Transmitting beam width (rad) 


io-< 


2 X 10-< 


10" 3 


Receiving aperture diameter (cm) 


o 


12.5 


15 


Laser amplifier 


None 


He-Xe 


None 


Detector 


Photo- 
multiplier 


In-Sb(77°K) 


Thermopile 


Measuring range (dB) 


65 


46 


50 
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the signal passes through a calcite polarizer, a quartz lens, and an iris, 
to a liquid-nitrogen-cooled indium-antimonide photovoltaic infrared 
detector whose output is fed into a synchronous detection system. The 
receiver acceptance angle, about 1 minute of arc, is determined by the 
dimensions of the amplifier tubes and the telescope. A signal-to-noise 
ratio of 45 dB is obtained on a clear day. The receiver sensitivity is 
limited by the laser amplifier noise. 

4.3 10.6fi Equipment 

A C0 2 laser with an output of 0.5 W and a thermopile at the focal 
plane of a 15 cm collecting mirror comprise the 10.6/x transmitting 
and receiving systems. The laser is operated with a continuous flow 
of gas. Therefore, virtually uninterrupted measurement has been pos- 
sible for a period of more than one year. 

4.4 3 mm Wavelength Radar 

The 3mm wavelength radar, based on a method used earlier for 
study of absorption by atmospheric gases at millimeter wavelengths, 20 
provides an independent measure of the liquid water content of fogs. 
The system uses a klystron, modulated linearly in frequency. Two 
precisely constructed trihedral corner reflectors are placed on the 
propagation path, one near the center of the path and the other at the 
optical transmitter terminal (see Figure 12). Measurement of the 
relative power reflected from the two reflectors results in attenuation 
data, accurate to about 0.2 dB per km, for that portion of the path 
between the reflectors. 

4.5 Rain Gauge System 

Three gauges of the tipping-bucket type which register a pulse for 
each 0.25 mm of rain, are located at the transmitting site, path mid- 
point and the receiving site (see Figure 12). Data from the three 
gauges are recorded at the site of the optical receivers. No quantita- 
tive measurements are made during snow. 

V. MEASURED RESULTS 

5.1 Fog 

In spite of available signal-to-noise ratios of 65 dB for 0.63/x, 45 dB 
for 3.5fi, and 50 dB for 10.6/x, the received signal at all three wave- 
lengths has often fallen below their noise levels during dense fog. Thus 
our ability to measure the wavelength dependence of attenuation by 
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fog is Limited to relatively light fogs with attenuations less than 25 
dB per km. Figure 13(a) shows that in a typical morning fog (Sep- 
tember 21, 1965) the attenuation (in decibels) of the 3.5/x wave is 
about half that of the 0.63/x wave. This result may be explained by 
fog droplets of about 1/x diameter. 

The measured data occasionally reveal a much smaller ratio as in 
Figure 13(b) where the 3.5/x attenuation (September 15, 1965) 
amounts to only one sixth of the 0.63/x attenuation; this result implies 
a fog with many droplets smaller than one micron. 

The measured attenuation of 10.6/x in a light fog is lower than at 
3.5/x, as shown in Figure 13(c). At the beginning of this record (1.30 
p.m. January 4, 1967), the signal at 0.63/x and 3.5/x are equal at a 
level of about —15 dB while the 10.6/x signal is at about —2 dB. 
This result can be explained, for example, by entering Figure 6(d) 
at a maximum drop radius of 1.3/*; there, one sees that the attenua- 
tions at 0.63/x and 3.5/x are equal and exceed the 10.6/x attenuation by 
a factor of eight, roughly the observed value. But at 3 p.m. in Figure 
13(c), the situation has changed significantly, and the ratios of the 
0.63/x, 3.5/x, and 10.6/x attenuations are now about 15:6:1. However, if 
we enter Figure 6(d) at a maximum drop radius of 0.7/x, the attenua- 
tion ratios are found to be very close to the observed values; thus it 
is found that small changes in the drop sizes can drastically affect 
wavelength dependence. 

A percentage time distribution for the ratio of attenuation by light 
fog at 3.5/x and 0.63/x is plotted in Figure 14. The data, taken during 
seven separate fogs, were accumulated over sixteen hours during 
which the attenuation at either wavelength over the 2.6 km path was 
between 5 and 50 dB. It is seen that the median value of the ratio is 
0.62 and the attenuation at 3.5/x is less than at 0.63/x 85 per cent of 
the time. The 15 per cent of the time that the attenuation at 3.5/x ex- 
ceeds that at 0.63/x may be caused in part by equipment instability; 
but it can be real for certain fog parameters, such as, in Figure 6(d), 
a m = 2 microns. 

The data of Figures 4 and 6 predict that attenuation of 10.6/x is 
insensitive to the drop-size distribution in light fog where drop radii 
are about 1/x. Therefore, one expects close correlation between the 
liquid-water density of a light fog and the measured 10.6/x attenua- 
tion. 

During measurement of attenuation of the laser beams by a morn- 
ing fog May 20, 1966 (Figure 15), the liquid water density of the 
fog was determined by measuring the attenuation of 3 mm waves. At 
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Fig. 13 — Measurement of 2.6 km transmission loss in light fog. dB-signal 
level in clear weather. 
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Fig. 14 — A distribution for the ratio of attenuation by light fog at 3.5ti and 
0.63m. 

3 mm, where the wavelength is certainly much greater than the fog- 
drop radius, the liquid water content is given by M = 0.25y gm/ni 8 , 18 
y being the 3 mm attenuation coefficient in dB/km. 

As shown in Figure 15, at 9:30 a.m. the laser beams could not be 
detected and the 3 mm data indicate a fog density of 0.3 gm/m 3 . As 
the fog gradually lifted, the optical signals emerge in sequence above 




Fig. 15 — Correlation between attenuation of laser beams and liquid water 
content measured by a 3 mm wavelength radar in fog. 
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the noise level but because dense haze lingered throughout the day, 
the 0.63/t wave remained below the clear weather level. There is good 
correlation between the attenuation of the laser beams and the liquid 
water content of the fog obtained by the 3 mm wavelength measure- 
ment. The zero reference of the liquid water density is taken at the 
vanishing of the attenuation of the 10.6/i. laser beam. However, precise 
agreement between the measured attenuation of 10.6/* and that calcu- 
lated from the liquid water content obtained at 3 mm is difficult to 
establish because of the very small attenuation at 3 mm for light fogs. 
It should also be kept in mind that the fog density measured by the 
3 mm radar is an average over only about half of the propagation 
path (see Figure 12). 

5.2 Rain 

Transmission through some light rains is difficult to interpret because 
the rain is mixed with fog ; however, summer showers are very suitable 
for evaluating attenuation by rain. Figure 16 shows typical measured 
data during a summer shower; attenuation of 3.5ti is consistently 
about 20 per cent greater than that of 0.63/x. The attenuation is well 
correlated with the rain rate averaged over the path. The average 
path rain rate is taken to be the arithmetic average of the readings 
of the three rain gauges described in Section 4.5, with that of the 
midpoint gauge carrying a double weight. 

Owing to the variation in raindrop terminal velocity and drop size 
distribution, the relationship between the attenuation coefficient and 
the rain rate can be hardly expected to be linear. However, the pre- 
dicted attenuation coefficient for the rain rate between 12.5 and 100 
mm per hour is almost linear as shown in Figure 8. Therefore, it is 
justifiable to determine a least square line for measured data of at- 
tenuation by rain above a rate of 12.5 mm per hour. Figure 17 shows 
a least square line 



•©-o^J + mi 



for the measured 0.63/x data of Figure 16 and two other showers. 

The calculated curves for 0.63/x attenuation by rain with and with- 
out forward scattering correction (from Figure 8) have been included 
in Figure 17 for comparison. The good agreement between the meas- 
ured line and the curve calculated with forward scattering may appear 
somewhat fortuitous in view of the rms deviation 2.8 dB/km. How- 
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Fig. 16 — Measurement of 2.6 km transmission loss in rain showers on Septem- 
ber 24, 1965. dB-signal level in clear weather. 



ever, only peak attenuations of summer showers which are relatively 
free of foggy conditions have been taken as measured data in Fig- 
ure 17. 

The relatively large rms deviation is explained as follows: the sam- 
pling of a spatially nonuniform shower by three rain gauges does not 
provide a satisfactoiy average rain rate over the 2.6 km path; the 
attenuation depends upon drop size distribution which may vary from 
one shower to another; since optical attenuation by fog is orders of 
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magnitude higher than that by rain, any slight mist or haze within 
the rain can appreciably change the measured values. 

Simultaneous measurements at all three wavelengths were carried 
out in an unusually heavy rainstorm of September 21, 1966, during 
which more than 125 mm of rain accumulated. A sample of the 
measured data from this storm is given in Figure 18. 




10 20 30 40 50 60 70 80 90 100 

p- AVERAGE PATH RAIN RATE IN mm/HR 

Fig. 17 — Measured rain attenuation coefficients of 0.63^ for a 2.6 km path; 
standard deviation is ±2.8 dB per km at a rain rate of 50 mm/hr. 



Points corresponding to the peak attenuations of the separate 
showers in the storm are plotted in Figure 19 and least-square lines 
for each wavelength have been drawn for rain rates between 10 and 
60 mm per hour. The attenuation by rain at 10.6/* and 3.5/i. is con- 
sistently greater than that at 0.63/*. At 0.63/* the agreement between 
calculated and measured values is again rather good. Notice that the 
measured points at 95 mm per hour fall below the predicted curve.* 
When we compare Figure 19 with Figure 13(c), say, it is obvious 
that the wavelength dependence in propagation through rain is in- 
verse to that through fog. 



* The 0.63/4 data in Figures 18 and 19 were not taken with the same receiving 
telescope as those of Figure 16 (5 cm diameter, see Table 1). Rather, a near-field 
cassegrain telescope with a collecting aperture of 600 cm was used. This change 
was made to determine whether or not a large receiving aperture would intercept 
more of the energy forward scattered by the rain (see Figures 10 and 21a) and 
thereby further decrease the attenuation. However, the attenuations (relative 
to the clear-day value) measured with the large aperture did not turn out to 
be significantly less than those measured with the small one. 
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Another sample of measured data is shown in Figure 20 ;f it is of 
particular interest because the attenuation at 0.63/* is much lower 
(about 18 dB) than in previous measurements and calculation for the 
same rain rates. The attenuation at 3.5/x and 10.6/* is only slightly 
below the prediction. Examination of the raw data shows that the 
rain distribution along the path was characterized by higher rates 




11 = 30 



Fig. 18 — Measurement of 2.6 km transmission loss in a rain storm. dB-signal 
level in clear weather. 



near the receiving site. This observation suggests that this very low 
attenuation at 0.63/x is caused by a strong forward scattering con- 
tribution which is most effective when the scattering particles are 
near the receiver. This result is also confirmed by the measured data 
of the rain showers in Figure 16. There the relatively low 0.63/x at- 
tenuation coefficient measured in the second shower is also associated 
with a rain distribution with high concentration near the receiving 
site. 



t The data in Figure 20 were taken with a 0.63/x transmitter beamwidth of 
5 X 10*' radians. 
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Fig. 19 — Measured attenuation coefficients of a rain storm for a 2.6 km path. 



5.3 Snow 

Figure 21 shows measurements of transmission loss through a light 
snow shower. The accumulation never exceeded about 12.5 mm per 
hour. Snows heavier than this have caused the signals at all three 
wavelengths to fall below the noise levels. The 10.6/a wave often ap- 
pears to suffer greater attenuation in snow as illustrated in Figure 21 ; 
however, no definite wavelength dependence is consistently observed. 

5.4 A Cumulative Distribution of Attenuation at 10.6 n 

In Section 4.3 we mentioned that the 10.6/* equipment operated 
continuously. Signal levels for the period June 1966 to May 1967 
have therefore been analyzed and reduced to a cumulative distribution 
of attenuation as shown by the central curve in Figure 22. For ex- 
ample, one sees that attenuation on the 2.6 km path, resulting from 
all causes, exceeded 30 dB (11.5 dB/km) three per cent of the time. 
Also shown in the figure are distributions for the months with the 
highest (March) and lowest (July) outage times. Clearly, during 
July, a system operating at 10.6/i, would be quite reliable; however, 
during March, heavy attenuations, mainly caused by fog, occur all 
too often on a 2.6 km path. 
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5.5 Beam Broadening at 0.63p 

Typical measured results on the beam broadening effects of rain, 
snow, and fog are presented in Figures 23, 24, and 25, respectively. 
Curve 1 in each figure is a typical clear-weather beam (about 20- 
second half-power beamwidth) ; this broadening is essentially deter- 
mined by the inhomogeneities of the refractive index in the atmos- 
phere. Notice that the clear weather beam itself is subject to modest 
variation but this variation is generally much less than the beam 
broadening observed under conditions of precipitation. The beam 




3:50 
E.D.S.T. (PM) 



Fig. 20 — Measured attenuation over 2.6 km during a rain shower on July 11, 
1967. dB-signal level in clear weather. 
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Fig. 21 — Measurement of 2.6 km transmission loss during a light snow shower 
on January 26, 1966. dB-signal level in clear weather. 
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Fig. 22 — A cumulative distribution of attenuation at 10.6/a over a 2.6 km 
path. O = Average attenuation from June 1966 to May 1967. A = March 1967, 
highest month. □ = July 1966, lowest month. 
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Fig. 23 — Beam broadening caused by rain. 
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Fig. 24 — Beam broadening caused by snow. 
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Fig. 25 — Beam broadening caused by fog. 

broadening is measured by moving the receiver on a horizontal track 
transverse to the axis of the incident beam. 

The other measured curves show beam broadening resulting from 
precipitation and the corresponding attenuations along with descrip- 
tions of precipitation; the curves are normalized at 6 = 0, the direc- 
tion of the beam axis. In general, the highest attenuations correspond 
to the greatest beam broadening. In Figure 24, we see that the onset 
of beam broadening caused by snow occurs when the attenuation is 
fairly high. Broadening effects in the case of rain are observable at 
much lower attenuations (Figure 23). This is because raindrops pro- 
duce more forward scattering than snow flakes which suffer from poor 
geometry and random orientations. Thus, the X points show the same 
amount of beam broadening in Figures 23 and 24; but the correspond- 
ing attenuations are 5.4 dB/km for the rain and 15.6 db/km for the 
snow. 

Although the beam broadening, caused by fog also is readily meas- 
urable, the fog drop, being comparable in size with wavelength, is even 
less efficient in forward scattering than the snow flake. The Y points 
in Figures 24 and 25 show that the attenuation by fog is about twice 
that of snow for the same degree of beam broadening. 

The measured beam broadening due to rain exceeds the prediction 
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of a single-scattering gaussian model as discussed in Section 3.2. 
Quantitative correlations are difficult because of the time variability 
and inhomogeneity of rainfall along the path. 

VI. DISCUSSIONS 

A simple picture of the vast difference in attenuation of micron 
waves by fog and rain has been substantiated by theory and experi- 
ment. The relatively low attenuation by rain is, in part, the result of 
forward scattering. Owing to the difficulty of obtaining independent 
physical data on snow, quantitative correlation between physical 
measurements of precipitation and transmission has been made only 
for rain and fog. The predictability of attenuation of optical waves 
by rain is comparable to that for microwaves. 21 Indeed, visible wave- 
lengths are attenuated less than millimeter wavelengths. 

Attenuation by dense fogs often exceeds 20 dB/km. The wave- 
length dependence of attenuation by light fog, that is, progressively 
less attenuation with increasing wavelength, has been observed as 
expected, but the dependence is a strong function of the drop-size 
distribution. Similar results for incoherent light have been reported 
in the literature. 1!5, 22 Although the predicted wavelength dependence 
for heavy fogs which consist of drops with average radius of 5/t, is 
weaker than for light fogs, the attenuation at 10.6/a would still be 
significantly less than at visible and shorter infrared wavelengths 
(see Figure 4). However, a typical dense fog of 0.1 g/m 3 liquid 
water density is calculated to attenuate 10.6/* by at least 40 dB/km. 
Attenuation by fog at 10.6/x is characterized by strong absorption 
and relative insensitivity to drop-size distribution. 

The measured attenuation by rain, about 0.2 (dB/km) /(mm/hr) 
at 0.63/x, is consistently less (about 20 to 40 per cent) than at 3.5/x 
and 10.6/i. Rain drops, which are absorbing spheres at 3.5/i and 
10.6/t, offer less scattering than at visible wavelengths and also 
have a less directional scattering pattern owing to the longer wave- 
lengths. 

The scattering patterns of fog drops have little directivity com- 
pared with the laser beams themselves therefore the forward scat- 
tering correction for the fog attenuation of narrow laser beams is 
small. If the beamwidth of the incident laser beam is comparable 
with the width of the scattering patterns of precipitation particles, the 
forward scattering will reduce the scattering loss significantly which, 
indeed, is mentioned in Section 5.2 for the propagation of 0.63/4 through 
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rain. Naturally the absorptive loss is not subject to any reduction. 
But the scattering loss of a laser beam propagating through fog can- 
not increase indefinitely as the fog density increases * In particular, 
for a visible source radiating into a homogeneously distributed dense 
fog, the energy must ultimately diffuse outward in all directions, thus 
the free space loss of an isotropic point source appears to be the 
limit of the scattering loss. 

When the scattering by precipitation particles is added to the at- 
tenuated original beam, a broadened beam is observed at the receiv- 
ing plane; this has been illustrated by a single scattering gaussian 
model in Section 3.2 and by measured results in Section 5.5. 

The apparent attenuation on the beam axis is reduced by the scat- 
tered component. For heavy rain, the received signal is dominated 
by the scattering as illustrated in Figure 10(b) and the reduction in 
attenuation becomes significant. When the scattering correction is 
small, estimates by the single scattering model (Section 3.2) and the 
multiple scattering model (Section 3.1) are practically the same. For 
relatively large corrections the multiple scattering model is expected 
to give a better estimate. 

The precise prediction of attenuation of laser beams by precipita- 
tion in the atmosphere is hampered by the uncertainty of particle 
size distributions. A better knowledge of the drop size distribution 
is certainly desirable to provide a firm ground for the prediction of 
attenuation, but even a measured drop-size distribution might lead to 
ambiguous interpretation when the variability and inhomogeneity of 
the natural fog or rain is unaccounted for. Summations over unreliable 
drop size distributions introduce complicated-looking problems but 
provide little more information than a simple relation such as equa- 
tion (10). 

Although it is convenient to distinguish between rain and fog by a 
classification of drop size, one should be aware of the frequent occur- 
rence of the intermediate states in the form of a drizzle or a mixture of 
rain and fog. The attenuation under these conditions has been observed 
to exhibit transition characteristics; however, little can be said beyond 
this trivial observation. 

The prediction of attenuation by snow has not been touched upon 
here, but the experimental data indicate that it lies intermediate be- 



*A common experience is to observe a heavy fog during daylight hours; al- 
though the sun is completely obliterated, the observer still has considerable 
scattered light in his environs. 
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tween that of rain and fog for a given amount of equivalent liquid 
water. 

In order to evaluate relative merits of systems using different wave- 
lengths, it is interesting to compare the attenuation by precipitation 
at optical wavelengths with that in other parts of the electromagnetic 
spectrum. Figure 26 presents the approximate attenuation coefficients 
for a typical dense fog and a representative shower in the frequency 
range from 10 10 to 10 15 cps. 18 - 21 The absorption for a water layer 
containing the same amount of liquid water as the rain is also shown 
in Figure 26 for reference. The attenuation by rain decreases slowly 
from millimeter to visible wavelengths. However, the attenuation by 
fog is about two orders of magnitude less than that of rain at upper 
microwave frequencies, while the optical attenuation by fog is about 
two orders of magnitude greater than that of rain. Notice that, even 
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at, lO/i, where attenuation by fog passes through a weak minimum 
(see Figure 26) , the attenuation is still rather large as far as a prac- 
tical communication system is concerned. 

ACKNOWLEDGMENT 

The authors are indebted to R. W. Wilson, A. A. Penzias, and 
A. W. Norris for the experimental data at 10.6/x and the distribution 
in Figure 22. Technical assistance was given by R. A. Desmond and 
C. A. Davison and R. A. Semplak; the computations were made by 
Mrs. E. Kerschbaumer. 

REFERENCES 

1. Chu, T. S. and Hogg, D. C, "The Attenuation of 3.392^ He-Ne Laser Radia- 

tion by Methane in the Atmosphere," B.S.T.J., 45 (February 1966), pp. 
301-306. , _ 

2. Middleton, W. E. K., Vision Through (he Atmosphere, University of Toronto 

Press, (1952) p. 43. ^ . , „ 

3 Hogg, D. C, "Scattering and Attenuation Due to Snow at Optical Wave- 
lengths," Nature, 203 (July 25, 1964), p. 396. p . 

4. Wilson, R. W. and Penzias, A. A., "Effect of Precipitation on Transmission 

Through the Atmosphere," Nature, 211 (September 3, 1966), p. 1081. 

5. H. C, Van de Hulst, Light Scattering by Small Particles, New York: John 

Wiley and Sons (1957), chapter 11. , 

6. Centeno, M., "The Refractive Index of Liquid Water in the Near Infrared 

Spectrum," J. Opt, Soc. Amer., 31 (March 1941), pp. 244-247. 

7. Bayly, J. G. and others "The Absorption Spectra of Liquid Phase H=0, 

HDO, and D 2 from 0.7/xm to 10/im," Infrared Phys., 3 (December 1963), 
pp. 211-223. . , , . 

8. Kislovskii, L. D., "Optical Characteristics of Water and Ice in the Infrared 

and Radio Wave Regions of the Spectrum," Opt. Spectroscopy, 7 (Septem- 
ber 1959), pp. 201-206. . 
9 Curcio J A. and Petty, C. C, "The Near Infrared Aborption Spectrum of 
Liquid Water," J. Opt. Soc. Amer.. 41 (May 1951), pp. 302-304 

10 Draegert, D. A., and others, "Far-Infrared Spectrum of Liquid Water, 

J. Opt. Soc. Amer., 56 (January 1966) , pp. 64-69. 

11 Eldridge, R. G., "A Few Fog Drop-Size Distributions," J. Meteorology, IS 

(October 1961), pp. 671-676. 

12. Diermendjian, D., "Scattering and Polarization Properties of Polydispersed 

Suspensions With Partial Absorption," Electromagnetic Scattering (Proc. 
Interdisciplinary Conf. Electromagnetic Scattering), ed. M. Kerker, New 
York: Pergamon Press, pp. 171-189, 1963. 

13. Erkovich, S. P., and others, "Influence of Fog on the Range of Ground Com- 

munication Using an Optical Carrier," Tele-communications and Radio 
Eng., (December 1965), pp. 12-16. trans, from Russian. 

14. Junge, C, "The Size Distribution and Aging of Natural Aerosols as Deter- 

mined From Electrical and Optical Data on the Atmosphere," J. Meteorol- 
ogy, 12 (February 1955), pp. 13-25. 
15 Kurnick, S. W., Zitter, R. N., and Williams, D. B. "Attenuation of Infrared 
Radiation by Fog," J. Opt. Soc. Amer., 60 (June 1960), pp. 578-583. 

16. Saunders, M. J., unpublished work. 

17. Laws, J. O. and Parsons, D. A., "The Relation of Rain Drop Size to Inten- 

sity," Trans. Amer. Geophys. Union, 2.', (1943), pp. 452-460. 



LASER BEAMS THROUGH PRECIPITATION 759 

18. Kerr, D. E., Propagation of Short Radio Wave, New York: McGraw-Hill 

Book Co., Inc., 1951, Chapter 8. 

19. Chu, T. S., "An Approximate Generalization of the Friis Transmission For- 

mula," Proc. IEEE, 53 (March 1965), pp. 296-297. 

20. Crawford, A. B. and Hogg, D. C, "Measurement of Atmospheric Attenuation 

at Millimeter Wavelengths," B.S.T.J., 35 (July 1956), pp. 907-916. 

21. Medhurst, R. G., "Rainfall Attenuation of Centimeter Waves: Comparison 

of Theory and Measurement," IEEE Trans. Antennas and Propagation, 
AP-13 (July 1965), pp. 550-564. 

22. Arnulf, A., and others, "Transmission by Haze and Fog in the Spectral 

Region 0.35 to 10 Microns," J. Opt. Soo. Amer., /,7 (June 1957), pp. 491-498. 



